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ABSTRACT:  Density  functional  theory  computations  were 
performed  to  investigate  the  stabilities  and  magnetic  and 
electronic  properties  of  VSe2  bulk,  few-layers,  monolayer,  and 
its  derived  nanoribbons  (NRs)  and  nanotubes  in  both  T  and  H 
phases.  All  these  materials  are  ferromagnetic,  but  exhibit 
versatile  electronic  properties.  The  VSe2  bulk  and  few-layers  in 
either  T  or  H  phase,  and  T  monolayer  are  metallic,  while  the  H 
monolayer  is  a  semiconductor.  For  nanoribbons,  the  zigzag 
NRs  in  both  T  and  H  phases  and  the  armchair  NRs  in  T  phase 
are  metallic,  while  the  armchair  NRs  in  H  phase  are  half- 
metallic.  The  edge  hydrogenation  can  retain  or  amend  the  electronic  property  of  the  pristine  NRs  depending  on  the  chirality  and 
phases.  Regardless  of  the  chirality,  nanotubes  in  T  phase  are  half-metallic,  while  those  in  H  phase  are  metals.  These  findings 
provide  a  simple  and  effective  route  to  tune  the  electronic  properties  of  VSe2  nanostructures  in  a  wide  range  and  also  facilitate  the 
design  of  VSe2-based  nanodevices. 


Bulk  Few-layer  Monolayer  Bare_NRs  II_NRs  NTs 


1.  INTRODUCTION 

Following  the  successful  exfoliation  of  graphene,1  many  two- 
dimensional  (2D)  layered  noncarbon  materials,2'3  such  as  BN 
and  MoS^  have  been  experimentally  realized  and  theoretically 
investigated.  These  2D  materials  are  bringing  revolutions  to 
numerous  advanced  applications  due  to  their  unique  and 
fascinating  physical  and  chemical  properties.  Moreover,  the 
properties  are  tunable  in  these  2D  materials  and  their  related 
one-dimensional  (ID)  nanoribbons,  for  example,  by  either  edge 
modification  or  surface  functionalizations,  and  by  rolling  them 
into  ID  nanotubes  with  different  chirality  and  diameters. 

Among  the  noncarbon  2D  materials,  transition-metal 
dichalcogenides  (TMDs)  have  been  extensively  studied.4  s 
The  layered  TMDs  consist  of  one  transition-metal  atom  between 
two  chalcogen  atoms  forming  covalently  bonded  sheets  with 
formula  MX,,  where  M  is  a  transition  metal  and  X  are  chalcogen 
atoms.  Adjacent  sheets  are  bonded  to  each  other  with  weak 
van  der  Waals  forces,  leading  to  highly  anisotropic  properties. 
The  anisotropy  of  the  materials  leads  to  interesting  properties 
connected  to  the  quasi  low  dimensionality,  and  therefore  TMDs 
serve  as  an  important  model  system  for  2D  phenomena. 
Typically,  the  properties  of  2D  nanomaterials  are  dependent  on 
the  thickness.  Taking  molybdenum  disulfide,  one  of  the  mostly 
studied  TMDs,  as  an  example,  the  MoS2  monolayer  exhibits  a 
~1.8  eV  direct  band  gap  while  a  bilayer  MoS2  exhibits  an  indirect 
band  gap  of  ~1.5  eV.9 

Vanadium  diselenide  (VSe2)  is  a  typical  TMD  material.10 
The  VSe2  monolayer  is  composed  of  metal  V  atoms  sand¬ 
wiched  between  two  Se  atom  (Figure  la),  and  the  Se— V— Se 
layers  are  stacked  in  the  (001)  direction  without  lateral 
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displacement  forming  trigonal  ( T )  phase  polytype  crystal 
(Figure  SI  a)  in  the  bulk.  Though  holding  very  similar  sandwich 
structure  with  MoS2  (Figure  lb),  VSe2  exists  only  in  T  phase 
(space  group:  P63mc)  in  nature,  while  MoS2  adopts  the 
hexagonal  (H)  phase  (space  group:  P63/mmc )  in  bulk  with  the 
Mo  atoms  in  octahedral  coordination  with  S  atoms. 

For  VSe2  bulk  crystal,  extensive  studies  have  been  performed 
to  investigate  its  charge  density  wave  (CDW)  transitions11  18 
and  intercalation  properties,19  “4  due  to  the  interaction  between 
the  electron  and  phonon  subsystem,25  and  the  ability  to  form 
intercalation  complexes  with  foreign  atoms  or  molecules 
incorporated  between  the  layers. 

Meanwhile,  the  VSe2  thin  films  have  been  synthesized  by 
different  approaches,  such  as  chemical  vapor  transport,26 
chemical  vapor  deposition,27'28  and  electron  beam  evaporation.29 
Quite  recently,  Xie’s  group  successfully  synthesized  2D  VSe2 
ultrathin  (four  to  eight  layers)  nanosheets  by  liquid  exfoliation  of 
a  bulk  VSe2  crystal  in  formamide  solvent,30  and  they  found  that 
such  nanostructures  retain  the  electronic  and  magnetic  characters 
like  its  3D  bulk  phase,  and  their  simulations  on  T-phase  bulk  and 
monolayer  reproduced  well  the  experimentally  observed  metallic 
properties. 

On  the  other  hand,  the  H-type  single-layer  components  were 
found  in  the  R  polytype  VSe2  bulk.  Spiecker  et  al.26  observed 
the  T  — *  R  polytype  transformation  in  the  thin  VSe2  surface 
layer,  which  was  induced  by  Cu  deposition.  Such  transformations 


Received:  July  16,  2014 
Revised:  August  12,  2014 

Published:  August  13,  2014 

dx.doi.org/10.1021/jp507093t  I J.  Phys.  Chem.  C  2014,  118,  21264-21274 


21264 


Figure  1.  Magnetism  distributions  (top  and  side  views)  of  2D  single-layer  VSe2  geometries  (3X3  supercell)  in  the  T  (a)  and  H  phases  (b),  as  well 
as  their  corresponding  hand  structures  (c,  d).  The  gray  and  khaki  atoms  represent  V  and  Se  atoms,  respectively.  The  isosurface  value  was  set  to  he 
0.008  e/A3.  The  spin-polarized  hand  structures  with  both  spin-up  and  spin-down  channels  are  presented  in  black  and  red,  respectively.  The  Fermi 
energy  is  denoted  by  a  green  dashed  line. 


were  also  found  when  elecrochemically  intercalating  Cu  into 
bulk  VSe231  or  intercalating  IC  into  VSe2  in  vacuum.  2  In  each 
single-layer  component  of  the  R  polytype  (see  Figure  Sid),  the 
V  atoms  are  in  trigonal  prismatic  coordination  with  Se  atoms, 
representing  a  H-type  monolayer  (Figure  lb).  Note  that  the 
only  difference  between  H-MoS2  (Figure  Sib)  and  R-V Se2 
(Figure  Sid)  bulk  is  the  stacking  patterns  of  the  sandwich 
layers:  the  adjacent  layers  were  stacked  by  tilting  the  crystal 
lattice  by  an  angle  (0)  of  ~17°  and  ~20°  from  the  c-axis 
orientation,  respectively. 

Exploring  the  VSe2  nanostructures  is  helpful  to  gain  insights 
into  the  properties  of  VSe2  at  the  atomic  scale.  So  far,  single¬ 
layer  VSe2  sheet  has  not  been  experimentally  realized;  a 
systematic  theoretical  investigation  on  its  structure  and 
properties  can  not  only  enhance  our  understanding  to  their 
intrinsic  characteristics  but  also  facilitate  their  applications. 
Despite  some  theoretical  studies  on  the  electronic  structures 
of  VSe2  bulk33,34  and  T-type  monolayer,30’35  many  important 
questions  are  pending:  What  is  the  stability  of  VSe2  bulk  in 
different  phases?  How  do  the  electronic  and  magnetic  properties 
vary  with  the  bulk-phase  structures?  What  is  the  stability  of 
T-  and  H-V Se2  monolayer  and  its  ID  derivatives?  How  do  the 
electronic  and  magnetic  properties  of  the  T-  and  H-V Se2  few- 
layers,  monolayer,  nanoribbons,  and  nanotubes  differ  from  the 
bulk?  How  can  chemical  modifications  such  as  hydrogena¬ 
tion  modulate  the  electronic  properties  of  VSe2  nanoribbons? 


In  this  work,  we  aim  to  address  these  issues  by  performing 
density  functional  theory  (DFT)  computations. 

2.  MODELS  AND  COMPUTATIONAL  METHODS 

Bulk  VSe2  forms  layered  crystal  (space  group  p63mc)  with 
the  separate  layers  stacked  along  the  (001)  direction  (c-axis 
direction  in  Figure  Sla).  Figure  1  presents  the  2D  monolayer 
models  used  in  this  work.  The  T  (Figure  la)  and  H  (Figure  lb) 
structural  constructions  are  the  same  as  TiS236  and  MoS2,37 
respectively. 

The  structural  models  of  the  studied  VSe2  nanoribbons  were 
constructed  from  the  2D  monolayer  slab  with  the  cleavaged 
plane  along  the  (001)  direction  (referred  to  as  the  c-axis 
direction  in  the  bulk  phase,  see  Figure  Si).  As  in  graphene  and 
other  inorganic  nanoribbons,  depending  on  the  cutting  direction, 
two  types  of  VSe2  ID  nanoribbons,  namely,  a -N  with  armchair 
edges  and  z-N  with  zigzag  edges  (N  is  the  ribbon  width),  can  be 
obtained.  The  width  N  is  represented  by  the  number  of  V  atoms 
in  a  unit  cell,  and  can  be  either  odd  or  even.  As  an  example, 
Figure  2  illustrates  the  top  and  side  views  of  a-I3  and  z-S  ribbons 
in  both  T  and  H  phases. 

Periodic  DFT  computations  were  performed  by  using  the 
plane-wave  technique  implemented  in  the  Vienna  Ab  initio 
Simulation  Package.38-4  The  ion— electron  interaction  is 
described  with  the  projector  augmented  wave42,43  method.  The 
exchange-correlation  energy  is  described  by  the  functional  of 
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Figure  2.  Geometries  of  the  top  and  side  views  of  ID  single-layered  VSe2  nanoribbons  in  IT  (left)  and  2 H  (right)  phases  with  armchair  (a,  b)  and 
zigzag  (c,  d)  edges  containing  two  unit  cells  without  optimization  [the  width  Na  =  13  in  (a)  and  (c),  Nz  =  8  in  (b)  and  (d)].  The  gray  and  khaki 
atoms  represent  V  and  Se  atoms,  respectively. 


Perdew,  Burke  and  Ernzerhof  (PBE)44  based  on  generalized 
gradient  approximation  with  spin  polarization.  A  500  eV  cutoff 
for  the  plane-wave  basis  set  was  used  in  all  the  computations. 
ID  periodic  boundary  condition  was  applied  along  the  growth 
direction  of  the  ribbons  to  simulate  the  infinitely  long  ribbon 
systems,  and  the  supercell  was  large  enough  to  ensure  a  vacuum 
spacing  greater  than  10  A.  For  geometry  optimization,  the 
k-point  grids  of  8  X  8  X  6,  8  X  8  X  1,  and  1X1X8  were  used 
for  3D  bulk,  2D  sheets,  and  ID  ribbons,  respectively,  while  larger 
k-point  grids  (25  X  25  X  1  and  21  k  points  for  VSe2  nanosheet 
and  nanoribbons,  respectively;  the  k  points  for  nanoribbons  are 
with  a  line  path  in  the  Brillouin  zone)  were  used  for  band 
structure  computations.  Note  that  the  standard  PBE  functional 
can  not  well  describe  the  van  der  Waals  interactions  between 
individual  layers.  Thus,  we  adopted  the  PBE-D2  method45  with 
the  inclusion  of  dispersion  energy  correction  for  bulk  and  few- 
layer  VSe2. 

For  T  bulk,  the  optimized  in-plane  (a  and  b  axis)  and 
interlayer  (c  axis)  lattice  parameters  are  3.33  and  6.17  A, 
respectively,  which  are  very  close  to  the  experimental  values 
of  3.35  and  6.12  A.46  The  T  bulk,  with  an  interlayer  separation 
(d)  of  3.04  A,  is  a  magnetic  metal  (Table  SI,  Figure  Sla),  in 
agreement  with  experimental  observation.47  For  comparison, 
we  also  optimized  the  bulk  VSe2  with  three  H  phases  (HI,  H2, 
and  H3 )  with  PBE-D2  method,  where  trigonal  prismatic 
monolayer  were  stacked  by  tilting  the  crystal  lattice  through  an 
angle  (6)  of  ~17°,  ~16°,  and  ~20°  from  the  c-axis  orientation, 
respectively  (Table  SI,  Figure  Sib— d). 

To  evaluate  the  stabilities  of  VSe2  sheets  and  nanoribbons, 
we  calculated  the  binding  energy  per  atom  (Fb)  using  the 
following  definition:  £b  =  (nEv  +  mFSe  —  Ev$sl)/(n  +  m)>  where 
Ev>  ESe,  and  EvSe,  represent  the  total  energies  of  free  V  atom, 
Se  atom,  and  VSe2  per  unit  cell,  respectively,  n  ( m )  denotes  the 
number  of  V  (Se)  atoms  in  a  unit  cell.  According  to  this 


definition,  systems  with  larger  Eh  values  (positive)  are  energeti¬ 
cally  more  favorable. 

For  the  edge-hydrogenated  nanoribbons,  the  average  binding 
energy  per  H  atom  is  defined  as  Eb  =  ( nEH  +  ENR  —  EH.NR)/n, 
where  EH,  £NR/FH.NR)  and  n  are  the  energy  of  single  H  atom, 
total  energy  of  pristine/hydrogenated  VSe2  NR,  and  number  of 
H  atoms,  respectively.  On  the  basis  of  this  definition,  systems 
with  stronger  binding  strengths  have  larger  (positive)  Eh. 

3.  RESULTS  AND  DISCUSSION 

3.1. 3D  VSe2  Bulk  with  Different  Phases.  Besides  Tbulk, 
we  also  computed  three  H  phases,  namely,  HI,  H2,  and  H3, 
for  VSe2  bulk.  The  T  bulk  has  an  interlayer  separation  (d)  of 
3.04  A,  while  the  corresponding  separations  between  layers  in 
HI,  H2,  and  H3  phases  are  3.13,  3.20,  and  3.25  A  (Table  l), 
respectively.  These  four  considered  bulk  phases  have  the  same 
stability  as  indicated  by  the  same  binding  energy  per  atom 
(5.22  eV).  Meanwhile,  our  computations  show  that  these  four 
phases  of  VSe2  bulk  are  all  ferromagnetic  metals  (the  magnetic 
moments  are  0.54,  2.68,  1.82,  and  2.72  //B  per  unit  cell  for 
T,  HI,  H2,  and  H3  phases,  respectively;  see  Figure  Si). 

3.2.  Pristine  2D  Monolayer  and  Few-Layer  VSe2 
Sheets  with  Different  Phases.  Very  recently,  Xie’s  group 
successfully  obtained  ultrathin  VSe2  nanosheets  with  thick¬ 
nesses  ranging  from  2.28  to  4.65  nm  (about  four  to  eight  Se— 
V— Se  atomic  layers),30  and  their  band  structure  computations 
revealed  that  the  T-type  monolayer  is  magnetic  metallic.  Here, 
we  studied  VSe2  monolayer  of  both  T  and  H  phases  (Figure  l), 
and  few-layer  VSe2  nanosheets  in  T,  HI,  H2,  and  H3  phases 
(Figure  Si).  The  number  of  Se— V— Se  atomic  layers  for  the 
few-layer  nanosheets  was  chosen  as  2,  3,  4,  and  8. 

For  VSe2  monolayer,  the  Se— V  bonds  in  T  and  H  phases  are 
2.48  and  2.50  A,  respectively,  and  the  corresponding  Se---Se 
distances  are  3.69  and  3.20  A.  The  thickness  of  the  T  mono- 
layer  is  3.15  A,  slightly  larger  than  that  in  bulk  counterpart  (3.13  A). 
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Table  1.  Calculated  Binding  Energies  per  Atom  (Eb,  in  eV), 
Interlayer  Separations  ( d ,  in  A),  and  Magnetic  Moments 
(M,  in  /<B)  of  T-,  HI-,  H2-,  and  H3-V Se2  Monolayer,  Two- 
Layer,  Three-Layer,  Lour-Layer,  Eight-Layer,  and  Bulk 


monolayer 

two-layer 

three-layer  four-layer 

eight-layer 

bulk 

T-V Se2 

Eh 

4.95 

5.18 

5.19 

5.20 

5.21 

5.22 

A 

3.04 

3.04 

3.04 

3.04 

3.04 

M 

0.44 

1.01 

1.53 

2.07 

3.12 

0.54 

Hl-V  Se2 

4.97 

5.19 

5.20 

5.21 

5.21 

5.22 

A 

3.13 

3.13 

3.13 

3.13 

3.13 

M 

1.00 

1.89 

2.78 

3.67 

7.23 

2.68 

H2-V  Se2 

Eh 

4.97 

5.19 

5.20 

5.21 

5.21 

5.22 

A 

3.20 

3.20 

3.20 

3.20 

3.20 

M 

1.00 

1.90 

2.80 

3.71 

7.35 

1.82 

H3-V  Se, 

4.97 

5.18 

5.19 

5.20 

5.20 

5.22 

A 

3.26 

3.26 

3.26 

3.26 

3.25 

M 

1.00 

1.89 

2.79 

3.70 

7.31 

2.72 

The  binding  energies  per  atom  for  T-V Se2  monolayer  (Eb, 
4.95  eV)  is  a  little  weaker  than  that  of  bulk  phase  (5.22  eV), 
while  the  Eb  of  the  H  monolayer  (4.97  eV)  is  slightly  stronger 
than  that  of  the  T  monolayer.  This  relative  stability  trend  is 


similar  to  the  case  of  the  VS2  monolayer.48  The  electronic 
structure  of  the  T  monolayer  resembles  that  of  the  bulk  very 
closely:  a  ferromagnetic  metal  with  the  magnetic  moment  of 
0.44  fi B  per  unit  cell  (Figure  la),  in  line  with  Ma  et  al.’s  recent 
theoretical  results.35  In  comparison,  the  slightly  more  stable 
H  monolayer  is  ferromagnetic  semiconducting  (l.OO  fiB  per  unit 
cell)  with  a  narrow  indirect  band  gap  (~0.19  eV,  Figure  lb), 
which  is  consistent  with  Pan’s  DFT  study.49  The  magnetism  of 
T  and  H  monolayers  originates  from  the  V  atom,  as  indicated  by 
the  magnetism  distributions  (Figure  la  and  lb). 

We  then  extended  our  studies  to  the  stability  and  magnetic 
and  electronic  properties  of  the  2D  VSe2  few-layer  nanosheets. 
The  interlayer  separations  ( d )  of  the  T-type,  HI -type,  and  H2- 
type  two-,  three-,  four-,  and  eight-layer  sheets  are  the  same  as 
that  in  bulk  (3.04,  3.13,  and  3.20  A,  respectively),  while  the 
d  value  of  the  H3-phase  few-layer  sheets  (3.26  A)  are  slightly 
enlarged  with  respect  to  that  in  their  corresponding  bulk  phase 
(3.25  A).  As  expected,  the  binding  energy  per  atom  (Lb) 
increases  as  the  number  of  Se— V— Se  atomic  layers  increases 
(Table  l),  but  the  Eb  values  of  eight-layer  nanosheets,  the 
thickest  few-layer  we  studied,  are  still  slightly  smaller  than  those 
of  bulk.  Among  the  four  phases,  the  two-layer  (Eb  =  5.18  eV), 
three-layer  (Eb  =  5.19  eV),  and  four-layer  (Eb  =  5.20  eV)  sheets 
in  T  and  H3  phases  have  lower  stabilities,  and  the  eight-layer 
sheets  in  T,  HI,  and  H2  phases  possess  the  same  binding 
energy  per  atom  (5.21  eV),  and  are  only  slightly  more  favorable 


Figure  3.  Magnetism  distributions  (top  and  side  views)  of  pristine  a-7  nanoribbons  (l  X  2  supercell)  in  the  T  (a)  and  H  phases  (b),  and  pristine  z-5 
nanoribbons  (1X2  supercell)  in  the  T  (c)  and  H  phases  (d),  as  well  as  their  corresponding  band  structures.  The  gray  and  khaki  atoms  represent  V 
and  Se  atoms,  respectively.  The  isosurface  value  was  set  to  be  0.008  e/A3.  The  spin-polarized  band  structures  with  both  spin-up  and  spin-down 
channels  are  presented  in  black  and  red,  respectively.  The  Fermi  energy  is  denoted  by  a  green  dashed  line. 


21267 


dx.doi.org/10.1021/jp507093t I J.  Phys.  Chem.  C2014,  118,  21264-21274 


Article 


The  Journal  of  Physical  Chemistry  C 


Figure  4.  Magnetism  distributions  (top  and  side  views)  of  pristine  a-8  nanoribbons  (1x2  supercell)  in  the  T  (a)  and  H  phases  (b),  and  pristine  z-6 
nanoribbons  (1X2  supercell)  in  the  T  (c)  and  H  phases  (d),  as  well  as  their  corresponding  band  structures.  The  gray  and  khaki  atoms  represent  V 
and  Se  atoms,  respectively.  The  isosurface  value  was  set  to  be  0.008  e/A3.  The  spin-polarized  band  structures  with  both  spin-up  and  spin-down 
channels  are  presented  in  black  and  red,  respectively.  The  Fermi  energy  is  denoted  by  a  green  dashed  line. 


than  H3  eight-layer  (Eb  =  5.20  eV).  The  magnetic  moments  of 
few-layer  sheets  of  these  four  phases  all  increase  with  increasing 
the  number  of  Se— V— Se  atomic  layers  (Table  l).  Similar  to  the 
electronic  properties  of  their  bulk  phases,  the  T,  HI,  H3,  and 
H3  few-layer  sheets  all  maintain  the  metallic  characteristics 
(Figures  S2— S5). 

Considering  the  comparable  stabilities  and  the  same  metallic 
electronic  properties  of  few-layer  nanosheets  in  the  four  phases, 
it  is  hard  to  say  that  the  experimentally  obtained  ultrathin 
VSe2  nanosheets30  with  metallic  character  is  in  T,  HI,  H2,  or 
H3  phase,  and  it  is  very  promising  to  synthesize  the  VSe2 
monolayer  and  nanoribbons,  which  may  exist  in  both  T  and  H 
phases. 

3.3.  ID  Bare  VSe2  Nanoribbons  (NRs).  Two  types  of 
nanoribbons  with  either  armchair-  or  zigzag-terminated  edge 
can  be  obtained  depending  on  the  cutting  directions  (Figure  2). 
Flere,  we  considered  both  odd  and  even  numbers  for  the  width 
N  and  examined  Na  =  7,  8  and  Nz  =  5,  6  as  representatives. 

Geometries.  The  optimized  a-7  VSe2  NRs  of  T  and  H  phases 
have  widths  of  10.81  and  9.77  A,  respectively.  The  edge  atoms 
in  armchair  ribbons  are  strongly  distorted.  Especially  for  the 
a-7  T-NR  (Figure  3a),  the  edge  Se— V  armchair  chains  on  both 
sides  are  optimized  to  a  zigzag  chain,  where  the  edge  Se  atoms 
shift  outward  (inward)  in  the  surface  plane  (edge  plane). 
Consequently,  the  bond  distances  of  edge  Se  to  the  edge  V  and 
subedge  V  are  elongated  to  2.25  and  2.67  A,  respectively 
(in  comparison,  the  Se— V  bond  length  in  the  2D  T-VSe2 
monolayer  is  2.48  A).  For  the  a-7  H-NR  (Figure  3b),  the  edge 


Se— V  armchair  chains  on  both  sides  retain  the  armchair  shape, 
while  the  edge  Se  atoms  shift  outward  in  the  surface  plane  and 
the  edge  V  atoms  are  drawn  close  to  the  edge  Se,  leading  to 
the  bond  length  of  edge  Se  to  subedge  (edge)  V  outstretched 
(contracted)  to  2.57  (2.30)  A,  where  the  Se— V  distance  is 
2.50  A  in  the  H-V Se2  monolayer  slab. 

The  widths  of  the  optimized  z-5  VSe2  NRs  of  T  and  H 
phases  are  13.35  and  12.18  A,  respectively.  Both  z-5  T-  and 
H-NRs  keep  the  zigzag  edges  (Figure  3c, d).  On  the  left  side, 
the  edge  Se— V  bonds  are  shrunken  to  2.32  and  2.44  A  for 
T  and  H  phases,  respectively,  while  the  subedge  Se— V  bonds 
are  shortened  (elongated)  to  2.36  (2.66)  A  for  z-5  T(H)-NR; 
on  the  right  side,  the  edge  and  subedge  Se— V  bonds  of  z-5 
T-NR  are  reduced  to  2.24  and  2.44  A,  respectively,  while  the 
bond  distances  of  edge  V  to  inner  Se  (parallel  to  width  direction) 
are  outstretched  to  2.84  A  (Figure  3c);  in  contrast,  the  edge  and 
subedge  Se— V  bonds  of  z-5  H-NR  are  slightly  elongated  and 
shortened  to  2.51  and  2.39  A,  respectively  (Figure  3d). 

The  VSe2  nanoribbons  with  even  number  V  atoms  (a-8  and 
z-6)  for  both  phases  have  very  similar  edge  reconstructions  with 
the  odd  numbered  NRs  (a-7  and  z-5),  as  shown  in  Figure  4. 
The  optimized  a-8  and  z-6  VSe2  nanoribbons  of  T  (H)  phase 
have  the  widths  of  12.56  (12.27)  and  16.34  (15.05)  A, 
respectively. 

Stabilities.  To  examine  the  feasibility  to  experimentally 
obtain  the  VSe2  NRs,  we  calculated  the  binding  energy  per 
atom  (£b)  as  the  indicator  to  evaluate  their  thermodynamic 
stabilities. 
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(a)  (b) 


Figure  5.  Magnetism  distributions  (top  and  side  views)  of  edge-hydrogenated  a-7  nanoribbons  (1X2  supercell)  in  the  T  (a)  and  H  phases  (b),  and 
edge-hydrogenated  z-5  nanoribbons  (1X2  supercell)  in  the  T  (c)  and  H  phases  (d),  as  well  as  their  corresponding  band  structures.  The  gray  and 
khaki  atoms  represent  V  and  Se  atoms,  respectively.  The  isosurface  value  was  set  to  be  0.008  e/A3.  The  spin-polarized  band  structures  with  both 
spin-up  and  spin-down  channels  are  presented  in  black  and  red,  respectively.  The  Fermi  energy  is  denoted  by  a  green  dashed  line. 


For  the  T-V Se2  nanorbbions,  the  armchair  NRs  are  more 
stable,  the  a-7  and  a-8  NRs  have  the  same  binding  energy  per 
atom  (4.89  eV),  larger  than  that  of  the  zigzag  NRs  (4.75  and 
4.78  eV  for  z-5  and  z-6,  respectively).  In  contrast,  for  the  H-V Se2 
nanoribbons,  the  zigzag  NRs  are  thermodynamically  more 
preferred  than  the  armchair  NRs  (the  Eh  values  for  z-5,  z-6, 
a-7,  and  a-8  NRs  are  4.84,  4.86,  4.76,  and  4.82  eV,  respectively). 

We  further  examined  the  stabilities  of  a-13  and  z-8  as 
representatives  for  the  NRs  with  larger  width.  As  expected, 
the  binding  energy  increases  as  the  ribbon  width  increases. 
The  binding  energies  of  a-13  and  z-8  T  (H)-NRs  are  4.92  (4.88) 
and  4.82  (4.89)  eV  per  atom,  respectively,  whose  widths  and 
magnetic  moments  are  20.90  (20.61)  and  22.17  (20.87)  A  and 
8.28  (4.51)  and  5.61  (5.39)  pB  per  unit  cell,  respectively. 

For  the  short-width  NRs  we  considered  above,  with  the  same 
width,  energetically  the  zigzag  NRs  of  H  phase  are  slightly  more 
favorable  than  those  in  H  phase,  while  the  armchair  NRs  prefer 
the  T  phase.  However,  the  binding  energy  difference  between 
the  armchair  and  zigzag  NRs  in  different  phases  will  diminish 
with  increasing  the  ribbon  width.  Thus,  we  expect  that 
the  experimentally  synthesized  nanoribbons  with  wider  width 
have  comparable  Eb  values,  and  may  exist  in  both  T  and 
H  phases. 

Magnetic  and  Electronic  Properties.  Regardless  of  the 
ribbon  width,  both  armchair  and  zigzag  NRs  of  T  and  H  phases 
are  ferromagnetic.  The  magnetic  moments  of  T  (H)  a-7,  a-8, 
z-5,  and  z-6  nanoribbons  are  2.34  (2.8l),  0.05  (3.65),  2.71 
(2.07),  and  3.25  (3.46)  pB  per  unit  cell,  respectively.  The  same 


as  the  monolayer,  the  ferromagnetism  of  VSe2  NRs  arises  from 
the  3d  orbitals  of  the  V  atoms  (Figures  3  and  4). 

Both  ID  armchair  and  zigzag  T  nanoribbons  retain  the 
metallic  character  of  its  2D  monolayer  as  indicated  by  the 
bands  crossing  the  Fermi  level  in  both  spin-up  and  spin-down 
channels  (Figure  3a, c  and  Figure  4a, c).  However,  the  electronic 
properties  of  H-V Se2  NRs  rely  on  the  cutting  ways  and  differ 
from  that  of  2D  monolayer:  the  armchair  NRs  are  half-metallic 
independent  of  the  width  (Figure  3b  and  4b),  while  the  zigzag 
NRs  display  a  metallic  behavior  regardless  of  the  ribbon 
width  (Figure  3d  and  Figure  4d)  as  the  T-V Se2  NRs  with  zigzag 
edges. 

3.4.  Hydrogenated  VSe2  Nanoribbons.  Since  the  edge 
atoms  in  the  bare  armchair  and  zigzag  NRs  are  coordinately 
unsaturated,  it  is  necessary  to  explore  how  the  magnetic  and 
electronic  properties  of  VSe2  nanoribbons  rely  on  edge 
termination.  Hydrogenation  is  an  efficient  way  to  saturate  the 
dangling  bonds,  and  we  investigated  the  edge-hydrogenated  VSe2 
nanoribbons  to  study  the  underlying  transformations  of  their 
electronic  and  magnetic  properties,  which  is  indeed  an  efficient 
computational  tool  and  could  be  a  guideline  for  controlling 
experimental  conditions  to  reach  the  desired  nanodevices.  In 
this  regard,  a-7,  a-8,  z-5,  and  z-6  NRs  were  chosen  as  the 
representatives  to  investigate  the  edge  hydrogenation  effect. 

When  the  edge  dangling  bonds  are  saturated  by  H  atoms, 
both  armchair  and  zigzag  NRs  in  either  T  or  H  phase  become 
less  distorted,  as  illustrated  in  Figures  5  and  6.  The  average 
binding  energies  of  H  at  the  edge  site  for  T  ( H )  a-7,  a-8,  z-5, 
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(a)  (b) 


Figure  6.  Magnetism  distributions  (top  and  side  views)  of  edge-hydrogenated  a-8  nanoribbons  (1X2  supercell)  in  the  T  (a)  and  H  phases  (b),  and 
edge-hydrogenated  z-6  nanoribbons  (1X2  supercell)  in  the  T  (c)  and  H  phases  (d),  as  well  as  their  corresponding  band  structures.  The  gray  and 
khaki  atoms  represent  V  and  Se  atoms,  respectively.  The  isosurface  value  was  set  to  be  0.008  e/A3.  The  spin-polarized  band  structures  with  both 
spin-up  and  spin-down  channels  are  presented  in  black  and  red,  respectively.  The  Fermi  energy  is  denoted  by  a  green  dashed  line. 


and  z-6  ribbons  are  1.88  (2.18),  1.90  (2.06),  2.55  (2.24),  and 
2.55  (2.07)  eV,  respectively. 

The  hydrogenated  nanoribbons  keep  the  ferromagnetic 
behavior  of  the  bare  nanoribbons.  For  T-NRs,  the  edge- 
hydrogenated  a- 7,  z-5,  and  z-6  NRs  retain  the  metallic  character 
of  the  pristine  NRs  (Figure  5a, c  and  Figure  6c),  while  edge 
hydrogenation  converts  the  metallic  a-8  NR  to  half-metallic 
(Figure  6a).  For  H-NRs,  hydrogen  saturation  transforms  the 
original  metallic  z-5  and  z-6  NRs  to  half-metals  (Figure  5d  and 
Figure  6d),  while  converting  the  pristine  half-metallic  a-7  and 
metallic  a-8  NRs  into  semiconductors  (Figure  5c  and  6c).  The 
partial  charge  density  distributions  of  the  edge-hydrogenated 
a-7  and  a-8  H-NRs  are  very  similar  (Figure  7a— d).  The  valence 
band  maximum  (VBM)  of  the  spin-up  (spin-down)  channel  is 
largely  composed  of  the  edge  side  (edge  and  subedge)  V  atoms, 
and  the  conduction  band  minimum  (CBM)  is  distributed  onto 
the  inner  (edge)  V  atoms. 

Note  that  the  3d  transition-metal  chalcogenides  possess 
subtle  electron  correlation  effects,  particular  on  electronic  and 
magnetic  properties.  Therefore,  we  also  used  PBE+U  (U  =  3  eV) 
method  to  examine  selected  VSe2  bulk  (T  and  HI  phases), 
monolayer  ( T  and  H  phases),  bare  (z-5  T-NR  and  z-6  H-NR), 
and  hydrogenated  nanoribbons  (a-8  T-NR  and  a-7  H-NR).  The 
PBE+U  lattice  parameters  are  larger  than  the  PBE-D2  and 
experimental  data,  for  example,  the  optimized  in-plane  ( a  and  b 
axis)  and  interlayer  (c  axis)  lattice  parameters  of  T  bulk  are 
3.46  and  6.32  A,  respectively,  while  the  corresponding  PBE-D2 
(experimental)  values  are  3.33  (3.35)  and  6.17  (6.12)  A. 
PBE+U  gives  the  same  binding  energies  per  atom  for 


T-  and  HI  -bulk/monolayer,  as  PBE  does.  Compared  with 
PBE  results,  PBE+U  predicts  the  same  relative  stabilities  and 
electronic  properties  for  bare  and  edge-hydrogenated  nanorib¬ 
bons  (the  band  structures  from  PBE+U  calculations  are  shown  in 
Figure  S6),  but  much  larger  magnetic  moments,  for  example,  the 
magnetic  moments  of  T  bulk,  bare  z-5,  and  hydrogenated  a-8 
T-NRs  are  predicted  to  be  1.18,  6.68,  and  8.22  //B,  respectively, 
while  PBE  yields  the  corresponding  values  of  0.54,  2.71,  and 
3.84  /ab.  Since  PBE  gives  the  same  trend  as  PBE+U  for  the 
stabilities,  electronic  band  structures,  and  magnetic  properties, 
we  used  the  PBE  data  for  discussions  in  this  work. 

3.5.  VSe2  Nanotubes.  Since  folding  nanoribbons  can  lead 
to  nanotubes,  we  also  investigated  the  stabilities  and  magnetic/ 
electronic  properties  of  VSe2  nanotubes,  for  which  (10,0)  and 
(5,5)  nanotubes  were  chosen  to  represent  zigzag  and  armchair 
tubes. 

The  binding  energies  of  (10,0)  and  (5,5)  VSe2  nanotubes 
in  T  phase  (4.76  and  4.72  eV  per  atom,  respectively)  are 
stronger  than  those  in  H  phase  (4.71  and  4.62  eV  per  atom, 
respectively),  but  lower  than  those  of  nanoribbons  because  of 
the  highly  strained  structures  in  VSe2  nanotubes. 

Compared  to  the  ferromagnetic  and  metallic  T  armchair/ 
zigzag  nanoribbons,  the  T  (5,5)  nanotube  with  the  diameter 
of  13.45  A  is  ferromagnetic  (10.52  /(B  per  unit  cell)  and  half- 
metallic  (Figure  8a);  the  T-(l0,0)  tube  (Figure  8b),  whose 
diameter  is  14.93  A,  is  also  a  ferromagnetic  half-metal  (20.85  //B 
per  unit  cell). 

In  contrast  to  the  versatile  electronic  characteristics  of 
H  nanoribbons,  both  armchair  and  zigzag  H  VSe2  nanotubes 
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Figure  7.  Partial  charge  densities  of  the  top  valence  band  (left)  and  bottom  conduction  band  (right)  for  edge-hydrogenated  H-V Se2  a-7  nanoribbon 
of  spin-up  (a)  and  spin-down  (b)  channels,  and  a-8  nanoribbon  of  spin-up  (c)  and  spin-down  (d)  channels. 


are  ferromagnetic  metals  (Figure  8c, d).  The  diameters  and 
magnetic  moments  of  H-( 5,5)  and  H-(l0,0)  VSe2  tubes  are 
13.22  and  14.97  A  and  5.61  and  3.15  //B  per  unit  cell, 
respectively. 

The  strain  energies  (Estr)  of  (5,5)  and  (10,0)  nanotubes  are 
6.89  (10.65)  and  11.  89  (15.40)  eV  in  T  (H)  phase.  We  then 
deduced  the  moduli  of  VSe2  nanotubes  following  the  pheno¬ 
menological  model  based  on  the  principles  of  the  classical 
theory  of  elasticity,  Estr  =  irYLh^/UR,  where  Y,  L,  and  h  are 
Yang’s  modulus,  length  of  the  nanotube,  and  thickness  of  the 
monolayer,  respectively.  The  Yang’s  moduli  of  (5,5)  and  (10,0) 
nanotubes  in  H  phase  (2.37  and  2.41  TPa,  respectively)  are 
larger  than  those  in  T  phase  (1.64  and  1.83  TPa,  respectively). 
The  calculated  Yang’s  moduli  of  VSe2  nanotubes  are  larger 
than  the  reported  values  (~1  TPa)  of  single-walled  carbon 
nanotubes,5  indicating  the  stronger  elasticity  of  VSe2  nano¬ 
tubes.  This  phenomenon  can  be  explained  by  the  higher 
energetic  expense  for  rolling  a  three-atomic  layer  VSe2  to  a  tube 
than  that  for  rolling  the  carbon  monolayer. 

Considering  the  high  strain  in  the  narrow  (5,5)  and  (10,0) 
nanotubes,  we  further  explored  the  VSe2  (10,10)  and  (17, 0) 
nanotubes  as  representatives  with  larger  diameters  (~19  A) 
to  gain  a  reliable  understanding  of  the  magnetic  and  electronic 
properties.  The  strain  energies  are  dramatically  reduced  (as 
listed  in  Table  2),  and  the  magnetic  moments  are  significantly 
increased,  while  the  Yang’s  moduli  are  slightly  decreased 
(except  for  T-(l0,10)  nanotube,  whose  Y  value  is  larger  than 
that  of  T-( 5,5)  nanotube).  The  (10,10)  and  (17, 0)  nanotubes 
in  either  T  (ferromagnetic  half-metals)  or  H  (ferromagnetic 
metals)  phase  have  the  same  magnetic  and  electronic  behaviors 
with  (5,5)  and  (10,0)  NTs  (see  Figure  S7). 


Table  2.  Calculated  Diameters  (D,  in  A),  Binding  Energies 
per  Atom  (Eb,  in  eV),  Magnetic  Moments  (M,  in  flB),  Strain 
Energies  (Estr,  in  eV),  and  Yang’s  Moduli  (Y,  in  TPa)  of  the 
Examined  VSe2  Nanotubes 


D 

Fb 

M 

£,tr 

Y 

H 5,5) 

13.45 

4.72 

10.52 

6.89 

1.64 

T-(10,10) 

19.61 

4.87 

21.19 

5.02 

1.72 

T-(10,0) 

14.93 

4.76 

20.85 

11.89 

1.84 

T-(l7,0) 

19.69 

4.88 

35.37 

8.08 

1.64 

H-(  5,5) 

13.22 

4.62 

5.61 

10.65 

2.37 

H-(10,10) 

18.92 

4.86 

17.81 

6.56 

2.09 

H-(10,0) 

14.97 

4.71 

3.15 

15.40 

2.41 

H-(  17,0) 

19.59 

4.87 

29.60 

10.45 

2.14 

4.  CONCLUSION 

In  this  work,  we  presented  detailed  DFT  study  on  the 
geometric  structures,  stabilities,  and  magnetic  and  electronic 
properties  of  VSe2  3D  bulk,  2D  monolayer,  and  few-layer 
nanosheets  in  both  T  and  H  phases,  as  well  as  their  derived  ID 
nanoribbons  and  nanotubes.  These  materials  are  all  ferromag¬ 
netic,  but  have  rather  versatile  electronic  properties  as 
summarized  in  Figure  9. 

The  3D  VSe2  bulk  in  either  T  or  H  (regardless  of  the 
stacking  patterns)  phase  has  the  same  binding  energy  per  atom, 
and  exhibits  ferromagnetic  metallic  character.  VSe2  few-layers 
and  T  monolayer  retain  the  ferromagnetic  metallic  properties 
as  their  bulk  structures;  the  H  monolayer  is  slightly  more 
favorable  than  T  monolayer  energetically,  and  is  a  ferromag¬ 
netic  semiconductor. 
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Figure  8.  Top  and  side  views  of  (5,5)  and  (10,0)  nanotubes  in  the  T  (a,  b)  and  H  phases  (c,  d),  and  the  corresponding  band  structures.  The  gray 
and  khaki  atoms  represent  V  and  Se  atoms,  respectively.  The  spin-polarized  band  structures  with  both  spin-up  and  spin-down  channels  are  presented 
in  black  and  red,  respectively.  The  Fermi  energy  is  denoted  by  a  green  dashed  line. 


T  Phase 


Bulk 


Few-laver  Monolayer  Bare  NRs  H  NRs 


NTs 


HPha.su 


Figure  9.  Relative  stabilities  and  electronic  properties  of  VSe2  bulk,  few-layer,  monolayer,  bare  and  hydrogenated  nanoribbons,  and  nanotubes  in 
T  and  H  phases.  M,  HM,  and  S  denote  the  metallic,  half-metallic,  and  semiconducting  properties.  For  NRs  and  hydrogenated  NRs  (H-NRs),  and  the 
larger  Z/A,  means  higher  stability  for  zigzag  and  armchair  nanoribbons. 


For  the  short-width  bare  NRs,  with  the  same  width,  zigzag 
NRs  are  slightly  preferred  energetically  in  the  H  phase,  while 
armchair  NRs  are  more  favorable  in  the  T  phase;  the  zigzag 
NRs  in  both  T  and  H  phases  and  the  armchair  NRs  in  T  phase 
are  ferromagnetic  metallic,  while  the  armchair  NRs  in  H  phase 
are  ferromagnetic  and  half-metallic.  The  edge  hydrogenation 
keeps  the  ferromagnetic  property  of  the  bare  NRs,  but  retain  or 
amend  the  electronic  property  of  the  pristine  NRs  depending 
on  the  chirality  and  phases. 

With  the  same  diameter,  VSe2  nanotubes  in  T  phase  are 
energetically  more  favorable.  Regardless  of  the  chirality, 


nanotubes  in  T  phase  are  ferromagnetic  and  half-metallic, 
while  those  in  H  phase  are  both  ferromagnetic  metals.  The 
triple-atomic-layer  VSe2  nanotubes  possess  larger  Young 
modulus  values  than  those  of  single-walled  carbon  nanotubes. 

The  fast-growing  success  in  isolating  various  2D  monolayer 
crystals  makes  us  highly  optimistic  that  monolayer  VSe2  and  ID 
nanoribbons  can  be  realized  experimentally  in  the  near  future. 
Though  no  H-V Se2  bulk  was  found  in  experiments,  the 
comparable/more  stable  monolayer  and  ID  nanoribbons  may 
coexist  in  the  H  form  with  those  in  the  T  phase,  and  possess 
versatile  electronic  properties. 
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